use TableA4.dta,clear
tab province,gen(prov) 

//column 1-4   
dprobit happy_d  direct_election2002 nomination2002 age age_sq male  ///
ethnic  mood  prov*, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) replace coefastr se
 
dprobit happy_d direct_election2002 nomination2002  age age_sq male  ///
ethnic  mood married ccp cadre_now cadre_past edu_year   health l_houincpc net_wealth working_hour unemployed ///
prov*, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) append coefastr se

//the following two columns restrict the obs to respondents who are heads of household or heads' spouses  

dprobit happy_d direct_election2002 nomination2002  age age_sq male ///
ethnic ccp cadre_now cadre_past edu_year child health mood ///
l_houincpc net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s ///
sp_edu_new health_s unemployed_s  prov* if H2_101<=2, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) append coefastr se

dprobit happy_d direct_election2002 nomination2002 dist_town2002 minority incpc98 vpop2002  disaster2002 age age_sq male ///
ethnic ccp cadre_now cadre_past edu_year child health mood ///
l_houincpc net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s ///
sp_edu_new health_s unemployed_s muchabove_aveinc above_aveinc muchbelow_aveinc ///
below_aveinc better_living worse_living expbig_increase ///
expsmall_increase exp_decrease  prov* if H2_101<=2, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) append coefastr se

//column5-8
dprobit happy_dn  direct_election2002 nomination2002   age age_sq male  ///
ethnic  mood prov*, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3)append coefastr se  


dprobit happy_dn  direct_election2002 nomination2002   age age_sq male married  ///
ethnic ccp cadre_now cadre_past edu_year  health mood l_houincpc ///
net_wealth working_hour unemployed  prov*, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3)append coefastr se  
 
//the following two columns restrict the obs to respondents who are heads of household or heads' spouses  

dprobit happy_dn   direct_election2002 nomination2002   age age_sq male  ///
ethnic ccp cadre_now cadre_past edu_year child health mood l_houincpc ///
net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s sp_edu_new ///
health_s unemployed_s  prov* if H2_101<=2, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) append coefastr se 


dprobit happy_dn   direct_election2002 nomination2002 dist_town2002 minority incpc98 vpop2002  disaster2002 age age_sq male  ///
ethnic ccp cadre_now cadre_past edu_year child health mood l_houincpc ///
net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s sp_edu_new ///
health_s unemployed_s muchabove_aveinc above_aveinc muchbelow_aveinc below_aveinc ///
better_living worse_living expbig_increase expsmall_increase exp_decrease ///
prov* if H2_101<=2, cluster(Village_code)
outreg2 using tableA4.doc, nolabel bdec(3) tdec(3) append coefastr se  
